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(54) Particle measurement by acoustic speckle 

(57) The Invention provides a metiiod and an appa- 
ratus for measuring the size, concentration and size of 
distribution of particles in a fluid by determining the 
characteristics of an acoustic speckle signal of the par- 
ticles and relating these characteristics to the size and 
concentration characteristics of the particles. The 
method and the apparatus are especially useful for the 
measurement of the size distribution of small particles 
in optically opaque liquids. 



< 

(O 

o> 

CM 

o> 
o 



LU 



Printed by Xerox (UK) Business Seivices 
2.16.7 (HRS)/3.e 



EP 1 092 976 A2 



Description 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

[0001 ] The present Invention relates to the measurement of size, concentration and size distribution of particles in 
fluids, and more particularly to the measurement of size, concentration and size distribution of particles in fluids by 
using acoustic speckle. 

10 

Description pf the Rg!qt?c| Art 

[0002] The measurement of the size and concentration of small particles in fluids is important in many industrial 
and diagnostic applications. Control or monitoring of operations such as crystallization, filtration, combustion, prepare- 
rs tion of reaction feed streams, degree of reaction and the like, often depend on the ability to measure the size and con- 
centration of small particles that are components of gas or liquid process streams. 

[0003] A number of methods of measuring size and concentration of small particles are known in the art. Surveys 
are provided, for example, in Principles, Methods and Application of Particle Size Analysis^ Syvitski, J.. (Ed.), Cam- 
bridge University Press, London, (1991) and by Allen. T, in Particle Size Measurement, Chapman and Hall, London 
20 (1990). 

[0004] Various optical methods that use coherent light -- lasers, for example - as well as incoherent light have been 
reported to be useful for particle measurement in transparent or translucent fluids and in situations where particle con- 
centration is low to moderate. Optical methods are reported that apply almost all wavelengths of light including wave- 
lengths that are smaller^ about the same as, or larger than the size of the particles to be measured. See, e,g, , Etzler, et 
25 al, Part. Part. Syst Char,, 12 (1 995). 

[0005] Other well known particle measurement techniques include phase doppler anemometry, microscopy with 
optical scanning and image analysis, sieving, and particle counting by optical and electrical counters. Although all of 
these methods are useful in certain applications, each has limitations that make ft difficult or impossible to use in certain 
situations. 

30 [0006] Particle measurement by the use of ultrasonic energy has been shown to have useful advantages in certain 
applications - for example, in streams that are optically opaque, or where particle density is high. The attenuation of 
sound as a function of frequency (Acoustic Attenuation Spectroscopy) has been used in the Ultrasizer™ device of Mal- 
vem Instruments Ltd., Worchester, U.K. In this method, a transducer Inputs an acoustic signal having a certain fre- 
quency range into the fluid sample and a second transducer receives the attenuated portion of the original signal. 

35 Mathematical modeling permits the calculation of particle size and size distribution from the measured attenuation 
spectra. 

[0007] Ultrasonfc measurement of particles is also described by Behrman et ai, On-line ultrasonic particle monitor- 
ing of brewing operations, MBAA Tech. Quarterly, 24:72 - 76 (1987); Bouts et ai, An evaluation of new asphaltene 
inhibitors: Laboratory study and field testing, JPT, 782 - 787 (1 995); de Boer et ai, Screening of crude oils for asphalt 

40 precipitation: Theory, practice and the selection of inhibitors, SPE Production & Facilities, 55 - 61 (1995); Dickinson et 
ai, Ultrasonic investigation of the particle size dependence of crystallization in n-hexadecane in water emulsions, J, of 
Colloid and Interface Science, 142(1 ):A 03 - 11 0 (1 991 ); and Holmes et ai, A wide bandwidth study of ultrasound veloc- 
ity and attenuation in suspensions: Comparison of theory with experimental measurements, J. of Colloid and Interface 
Science, 755:261 - 268 (1993); among others. U.S. Patent Nos. 4,412,451 to Uusitalo etai, 4,509,360 to Erwin efa/., 

45 5,569,844 to Soweiby. and 4.706,509 to Riedel, also cover aspects of particle measurement by using ultrasonics. 
[0008] In U.S. Patent Application Serial No. 08/947,821 , several of the present inventors described a novel use of 
uftrasonics to measure the size and concentration of small particles in fluids. In that method, a focused acoustic signal 
was transmitted into a fluid. Acoustic energy that is scattered by particles in the focal region of the signal is sensed by 
the same or another transducer. The signal of scattered ultrasonic energy is transformed from a time fonnat to a fre- 

50 . quency format and particle size and size distribution can be calculated by analysis of the magnitude of the signal as a 
function of frequency. This method was reported to be useful for the real-time measurement of particles in opaque liq- 
uids, such as crude oil. 

[0009] A characteristic that has long been recognized in both laser and acoustic spectrometry is the phenomenon 
of "speckle". General discussions of ultrasonic speckle have been provided by Wagner etai, Statistics of speckle in 
55 ultrasonic B-scans, IEEE Transactions on Sonics and Ultrasonics, 30(3): 1 66 - 1 60 (1 983), and Abbott et ai, Acoustic 
speckle: Theory arui experimental analysis, Ultrasonic Imaging, 1 :303 - 324 (1 979). 

[0010] The presence of acoustic speckle has been traditionally regarded as an annoyance in applications such as 
ultrasonic medical imaging and efforts to minimize ultrasonic speckle have been reported by Wells et ai, Speckle in 
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ultrasonic imaging, Ultrasonics, 225 - 229 (Sept. 1981), and Bamber et aL, Ultrasonics, 41 - 44 (Jan. 1986), and 
speckle minimization lias been the subject of U.S. Patent Nos. 4,561 ,01 9 (RE35,1 48) to Lizzi ef a/., 4,771 ,470 to Geiser 
eta!., and 5,090,412 to Stiimazaki. 

[0011] More recently, however, there have been reports of the potential of deriving useful information from speckle 
5 patterns. See, e.g.. Hong, S-K, and J-B Han, Considerations on speckle pattern interferometry of ultrasonic speckles, 
Ultrasonics, 35:329 - 332 (1997), for a description of the potential for using ultrasonic speckle to produce information 
about defonnations and flaws in the internal parts of a material that can transmit ultrasonic waves. The same group 
report how known deformations produce changes In ultrasonic speckle patterns in Hong, S-K and Y G. Ohr, Ultrasonic 
speckle pattern correlation interferometry using a pulse-echo method, J. Phys. D: Appl. Rhys., 37:1392 - 1396 (1998), 
10 [001 2J The use of ultrasonic speckle signals to derive infomiation about particles In fluid systems has been 
described by Nakajima etal. (U.S. Patent No. 4,944,1 89) who report an ultrasonic speckle velocity measurement appa- 
ratus and method. A reported advantage of the method is that it can measure fluid velocity even at low velocities and 
the measurement is reported to be unaffected by the number or concentration of particles in the fluid. 
[001 3] Despite significant improvements in the methods and instruments available to measure the size and concen- 
ts tration of small particles In fluids, there remain certain applications where such measurement is difficult or impossible. 
For example, in fluids that are optically opaque - such as crude oil - or where particle concentration is high, and where 
measurements are required on a rapid basis - even on a real-time basis and without dilution of the fluid, conventional 
particle measurement systems are unsuitable. Accordingly, it would be useful to provide a method and an apparatus for 
measuring particle size, concentration and size distribution of small particles In fluids that would be suitable for use even 
20 in optically opaque fluids. It would also be useful to provide such a method and an apparatus that could provide such 
measurements on a rapid, or real-time basis. It would also be useful if such a method and apparatus could perform 
these measurements without diluting the fluid in which the particles are carried. 

Brief Summary pf the Invention 

25 

[0014] Briefly, therefore, the present invention is directed to a method for detenmining the size of particles in a fluid 
comprising the steps of obtaining an acoustic speckle signal of the particles in the fluid and deriving the size of the par- 
tteles in the fluid from the acoustic speckle signal. 

[0015] The present invention is also directed to a novel method for measuring the concentration of particles in a 
30 fluid comprising the steps of obtaining an acoustic speckle signal of the particles in the fluid and deriving the concen- 
tration of particles in the fluid from the acoustic speckle signal. 

[0016] The present invention is also directed to a novel method for determining the size distribution of particles in 
a fluid comprising the steps of obtaining an acoustic speckle signal of the particles in the fluid and deriving the size dis- 
tribution of the particles from the acoustic speckle signal. 
35 [001 7] The present invention is also directed to a novel method for determining the degree of particle agglomeration 
in a liquid containing particles that produce individual speckle events, the method comprising obtaining, for two or more 
individual speckle events, an individual acoustic speckle signal from the liquid and determining the degree of agglom- 
eration from the acoustic speckle signal. 

[0018] The present invention is also directed to a novel apparatus for determining the size distribution of particles 
40 in a fluid comprising a signal generator capable of generating an electrical signal; a transducer that can translate the 
electrical signal from the signal generator into an acoustic input signal that is introduced into the fluid; a sensor that can 
sense acoustic energy that is scattered by the particles In the fluid in response to their interaction with the acoustic input 
signal and provide an output electrical signal that is related to the sensed acoustic energy; an oscilloscope that is capa- 
ble of acquiring the output signal in an amplitude versus time domain and transfomriing the signal to a magnitude versus 
45 frequency domain; and a computer that is programmed with an algorithm that is capable of identifying an acoustic 
speckle signal in the frequency domain and relating the magnitude and/or duration of the acoustic speckle signal to the 
size distribution of the particles in the fluid, 

[0019] Among the several advantages found to be achieved by the present invention, therefore, may be noted the 
provision of a method and an apparatus for measuring particle size, concentration and size distribution that are suitable 
50 for use even in optically opaque fluids; the provision of such a method and an apparatus that could provide such meas- 
urements on a rapid, or real-time basis; and the provision of such a method and an apparatus that could perfomn such 
measurements without diluting the fluid in which the particles are carried. 

Brief Description of the Several Views of the Drawings 

55 

[0020] 

Figure 1 shows conceptual representations of an acoustic wave signal approaching particles in a fluid that are 
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arranged in (a) one fonm of an acoustic mirror, or (b) a Bragg lattice; 

Figure 2 shows a persistence display for over 1000 sweeps from an oscilloscope of an amplitude versus time plot 
. of acoustic energy scattered by particles in crude oil, showing (a) no speckle; (b) several speckle events within a 
focal region; and (c) speckle events within and outside of a focal region; 

Figure 3 is a block diagram of one embodiment of an apparatus suitable for practicing the present invention show- 
ing the arrangement of components; 

Figure 4 shows a representation of a cross-sectional view of a focused transducer of the present invention mounted 
in a pipe containing a flowing fluid and indicating the convergence of the acoustic signal into a focal region within 
the fluid; 

Rgure 5 is a plot of an acoustic speckle signal from crude oil after transformation by Fast Fourier Transform (FFT) 
into magnitude versus frequency domain, and of the same type of FFT data averaged over time, wherein an acous- 
tic speckle peak is apparent at a frequency of about 10-12 MHz. which is the frequency of the exciting transducer; 
Figure 6 is a series of plots representing an oscilloscope screen displaying a typical scattered acoustic energy sig- 
nal from crude oil during several steps of the method of the invention, where 6(a) shows the amplitude of sensed 
scattered acoustic energy in a time domain; 6(b) shows the amplitude of sensed scattered acoustic energy in a time 
domain and indicating the selection of that portion of the signal emanating from the focal region as a bolded portion 
of the signal; 6(c) shows the selected amplitude versus time signal of the signal of 6(b) after Fast Fourier transfor- 
mation into a frequency domain; and 6(d) is a display of the average of the Fast Fourier Transformed signal in a 
frequency domain over 900 sweeps with identification of a speckle peak as the broad peak on the left side of the 
display; 

Figure 7 is a representation of a plot of the frequency of occun-ence of speckle events produced by particles having 
an average size of a^vG versus the average size of particles (a avg) including one sub-range region; 
Figure 6 shows an oscilloscope display of an averaged transformed output signal from a crude oil sample that has 
a speckle peak at a .frequency of about 10-12 MHz superimposed with an averaged transformed output signal from 
a control fluid - n-heptane; 

Figure 9 shows a plot of speckle height that is derived from the difference between the averaged transformed out- 
put signals shown in Figure 8; 

Rgure 1 0 showis a conceptual view of an amplitude (A) versus time (6) plot of scattered acoustfc energy where the 
scan has been repeated at a selected repetition frequency over a period of time along the G' axis; and 
Figure 1 1 is a plot of the maximum height of the peak on each of the A vs. G plots of Rgure 10 plotted versus time 
(G'), and illustrates the duration of a speckle event. 

[0021] Con^sponding reference characters indicate corresponding parts throughout the several views of the draw- 
ings. 

Detailed Description of the Invention 

[0022] In accordance with the present Invention, it has been discovered that the size, concentration and size distri- 
bution of small particles in a fluid can be derived from a measurement of acoustic speckle signals produced by the par- 
ticles in response to input signals of acoustte energy. This surprising discovery is especially significant because 
acoustic speckle, particularly for energy of ultrasonic wavelengths, can be measured very quickly and can be transmit- 
ted and sensed even in optically opaque fluids. Moreover, the novel measurement technique of this invention can be 
used alone or can be combined with other particle measurement techniques to provide valuable knowledge about the 
size and concentration characteristics of both small and larger particles in a fluid on a real-time basis. The method of 
this invention can be used on fluid samples in a test cell or on fluids in process flow streams, in tanks or in reactors. In 
fact, it is believed that the method can be applied anywhere an ultrasonic transducer can be introduced into a fluid to 
be tested. Moreover, it is believed that the method can be applied to any fluid that will transmit ultrasonic waves and that 
it can be used to measure the size and concentration characteristics of any particles that scatter ultrasonic energy and 
that are sufficiently small that their motion within the fluid is controlled by Brownlan motion. 

[0023] Without wishing to be bound by this or any other theory, it is believed that the production of acoustic speckle 
signals by small particles in a fluid is related to their tendency to move randomly with Brownian motion. When a beam 
of monochromatic acoustic energy having a certain wavelength and a certain diameter interacts with the particles, the 
particles scatter some portion of the acoustic energy. From time to time, some of these particles, each having a diam- 
eter that Is less than the wavelength of the interacting beam of acoustic energy and less than the beam diameter, 
assemble stochastically into a common plane before moving out of formation again. Occasionally, this plane will lie 
orthogonal to the path of the acoustic beam and the backscattered signal from each particle will form a coherent sum 
at the transducer. When these criteria are satisfied, the particles are said to form an acoustic mirror, one form of which 
is illustrated in Figure 1 (a). A similar effect can also be produced by particles that are not in the same plane orthogonal 
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to the propagation direction of the ultrasound signal, as long as the particles lie in planes that are separated an integral 
number of wavelengths apart along the path of the sound wave. Such an arrangement is termed a "Bragg lattice", as 
illustrated in Figure 1 (b), which also results in a coherent sum of the backscattered signal from all the particles at the 
transducer. Particles in either of the arrangements described above will be referred to herein as being in an acoustic 
5 mln-or, or simply, a mirror. 

[0024J As used herein, particles are considered to be "small particles" if their motion in a static fluid is controlled 
' primarily by Brownian motion, as opposed to, for example, gravitational force. The diameter of such small particles is 
believed to be much less than the wavelength of the acoustic energy signal. For a signal having a wavelength that cor- 
responds to a frequency of 1 MHz, for example, this would mean particles having a diameter of from about 5 microme- 
10 ters down to a few nanometers. 

[0025] It would also be possible for larger particles to create coherent backscattered signals that sum at the trans- 
ducer as a result of mechanical stirring, convection currents, sedimentation or creaming. 

[0026] For the purposes of this invention, it will be assumed that all particles and agglomerates of particles are sub- 
stantially spherical particles having the same mass as the actual particle and the term "particle size", as that term is 
15 applied to single particles, is to be understood to mean a number that Is proportional to the diameter of the substantially 
spherical particle or the agglomerate. 

[0027] During the period that the scattering particles remain in an acoustic mirror, the approaching acoustic wave 
encounters each particle in the lattice at the same moment and acoustic energy is emitted by each scattering particle 
simultaneously and in phase. A portion of the acoustic energy that is scattered by these particles can be sensed and 

20 measured by a sensing transducer to produce an electrical signal that corresponds to the scattered acoustic energy 
The acoustic emissions from the min-or combine at the sensing transducer and produce a characteristic signal that is 
termed an "acoustic speckle" (which, hereinafter, may also be temned a "speckle", or "speckle signal", or "speckle 
evenr). The speckle can be displayed on an oscilloscope amplitude versus time plot of scattered ultrasonic energy as 
a peak having an amplitude significantly greater than that of the energy scattered by particles that are not in a mirror. 

25 Several acoustic speckle peaks can be seen in Figures 2(b) and 2(c) as dotted line excursions whose amplitude signif- 
icantly exceeds the amplitude of non-speckle scattering, 

[0028] If it is assumed that the formation and subsequent disintegration of the min-or Is due to the random motion 
of the particles, or aggregates of particles. In the fluid, then, by measuring certain properties of the speckle signal, infor- 
mation can be gathered regarding the mirror itself, and the particles contained in it. Although each mirror forms and dis- 

30 integrates due to the random motion of the scattering particles, the period of time for which the structure exists is 
dependent upon the speed with which the scattering particles move in the fluid. The velocity of small particles in a fluid 
due to Brownian motion depends upon the particle size and the viscosity of the fluid, in a fluid having a given viscosity, 
smaller, dispersed particles will move faster than their larger, or aggregated, counterparts and speckles from small par- 
ticles will appear and disappear more quickly than will speckles produced by large particles or aggregates. 

35 [0029] Surprisingly, it has been found that the duration of an acoustic speckle event can be used to derive the size 
of the particles in the minror that produced the speckle. The duration of the acoustic speckle event represents the life 
span of the speckle event in time. As used herein, the term "duration" is a number that is proportional to the life span 
overtime of a measured signal of a speckle event. Depending upon how the duration of a speckle event is measured, 
the measured value can be equal to the total life span of the individual signal of the speckle event, or it can be some 

40 particular fraction thereof. This concept can be described, for illustration purposes only, by reference to Figures 10 and 
1 1 . Rgure 1 0 shows a conceptual view of an amplitude (A) versus time (6) plot of scattered acoustic energy where the 
scan has been repeated at a selected repetition frequency over a period of time along the 9' axis. During the period of 
time that is shown, an acoustic speckle signal is sensed and the speckle signal is shown as the peak in the A vs. 6 plots. 
The maximum height of the peak on the A vs. 8 plots can be plotted versus time (6') as shown in Rgure 1 1 to Illustrate 

45 the duration of the speckle event. 

[0030] As described above, the duration of this speckle event can be derived from either a direct measurement of 
the life span on the persistence display of the oscilloscope, or correlation with frequency dependent scattering ampli- 
tude via scattering theory by using techniques analogous to those employed by dynamic light scattering instrumentation 
such as the Malvern Zetasizer, (Malvern Zetasizer Manual, Malvern Instruments, Malvern, U.K.). Alternatively, the dura- 

50 tion can be measured by the calculation of some parameter from the speckle signal peak that is a particular fraction of 
the life span -- such as by calculation of an autocorrelation function as described by Berry et al. , in Physical Chemistry, 
p.1106, Wiley Publ. New York (1980); and by Wellstead, P.E., I^ethods and Applications of Digital Spectml Techniques, 
Technical Report No. 0008/83, UMIST, Solartron Instruments, Farnborough, U.K.. Such an approach is illustrated In 
Figure 1 1 by the straight line passing through the maximum of the peak and through the curve at a selected later point 

55 The slope of the line can be used to derive the duration, since it can be related to the life span of the speckle event. 
[0031] It is believed that a particle size derived from the duration of a speckle signal by any of these methods rep- 
resents an average of the sizes of the particles that produced the speckle signal. When the term "particle size" is used 
herein, it is to be understood to include absolute particle size, as well as values that are related to absolute particle size 
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and from which the absolute particle size can be derived. 

[0032] It also has been found, surprisingly, that a measurement of the frequency of occun-ence of acoustic speckle 
events can be used to derive the concentration of particles in the fluid. It is believed that there is a direct con-elation 
between the concentration of the particles in the fluid and the number of speckle events that occur per unit time. The 
5 higher the particle concentration in the fluid, the more frequent the particles fomri a mirror and, thus, produce an acous- 
tic speckle event. 

[0033] Furthermore, it has been found that the size distribution of the particles in the fluid can be derived from a 
measurement of the number and duration of speckle events over a sufficient period of time. This derivation is possible 
because the average size of the particles producing each speckle event can be derived from the duration of the individ- 
10 ual speckle signal and, thus, the frequency of occun'ence of speckle events for an array of average particle sizes essen- 
tially represents the particle size distribution of the particles in the fluid. Since the measurement of speckle duration 
provides an average particle size, it is believed that the more speckle events that are included in the calculation, the 
more accurate will be the particle size distribution that is derived in this manner. 

[0034] Also, it has been found, surprisingly, that the peak height of an acoustic speckle signal is related to the total 
15 cross-section area of the particles that produced the speckle event. This provides the basis for an alternative method 
of deriving particle concentration. The method can be briefly described as follows. Since the average size of the parti- 
cles producing the speckle event can be derived from the duration of the speckle signal - as described above ~ the 
number of particles that produced the speckle can be derived by dividing the total cross section area of partrcles pro- 
ducing the speckle by some function of the average particle diameter. This number can be used to produce a particle 
20 concentration value when it is related to some characteristic volume, such as, for example, the volume of the focal 
region of the transducer, as will be described below. Once particle concentration is determined, the size distribution of 
the particles in the fluid can be derived by a method that is similar to that described above. 

Apparatus: 

25 

[0035] A system that may be used to practice the present invention is shown in Figure 3. In general, a signal gen- 
erator, or pulser 30, generates an electrical signal which is transmitted to a transducer 10 that transduces the electric 
signal into an acoustk: input signal that is introduced into the fluid having the particles to be measured. The same trans- 
ducer, or, optionally, a separate transducer, acts as a sensor 20, to sense a portion of the acoustic energy that is scat- 

30 tered by the particles in the fluid in response to their interaction with the acoustic input signal and provide an output 
electrical signal that con-esponds to the sensed acoustic energy An amplifier 40 amplifies the output signal and trans- 
mits the amplified signal to an oscilloscope 50 that is capable of acquiring the output signal in an amplitude versus time 
domain and transforming the signal to a magnitude versus frequency domain. If the transducer has a focal lens that 
focuses the input signal on a focal region in the fluid, the oscilloscope is preferably capable of selecting that part of the 

35 detected scattered signal that results from acoustic scattering by particles in the focal region (which step can be tenned 
"gating") and transfonning only the selected part of the output signal into a magnitude vs. frequency domain. This trans- 
formed output signal can be transmitted to a computer 60 that is programmed with an algorithm that is capable of relat* 
ing the magnitude and/or duration of the acoustic speckle signals to the size of the particles in the fluid. 
[0036] As would be readily recognized by one of ordinary skill in the art of ultrasonic measurement, the system 

40 described above could be easily modified while still carrying out the same functions. For example, the electric signal to 
be transmitted to the transducer 10 could be generated by a combined pulser/amplifler, as well as by the pulser 30. 
Alternatively, the oscilloscope 50 could provide amplification of the sensed signal and replace the separate amplifier 40. 
Likewise, if desired, the computer 60 could accomplish the calculation functions ascribed above to the oscilloscope 50. 
or vice versa. The acoustic transducer 10 may be mounted in a test cell, or may be mounted directly in a pipeline or a 

45 tank. Other systems for the ultrasonic measurement of particles in liquids are generally described in U.S. Pat. No. 
4.412,451 to Uusitalo ef a/., U.S. Pat. No. 4,509,360 to Erwin ef a/., U.S. Pat No. 4.706,509 to Riebel, and U.S. Pat. No. 
5,121.629 to Alba. 

[0037] The various parts of a system to practice the present invention and their operation are described as follows. 
[0038] The transducer 10 is an ultrasonic transducer that introduces a signal of acoustic energy into the fluid The 

50 transducer includes, in general, a piezoelectric element capable of transforming electrical signals into physical pulses. 
If the element Is in contact with a fluid, such physical pulses are transfen-ed to the fluid and initiate waves having a fre- 
quency that is conti-olled by tiie frequency of the electronic signal. Preferably, the transducer 10 is a focused transducer 
that contains a lens to focus the signal as illustrated by the dashed lines of Figure 4. The focal length of the transducer 
is the distance from the end of the transducer to the point where the wave patterns converge. Focal length may be esti- 

55 mated as: 

Focal length = R (C ^)/C , - C ^) 



6 



EP 1 092 976 A2 



where: R is the radius of the lens of the transducer; is the speed of sound in the lens; and is the speed of sound 
in the probed medium. 

[0039] While focal length is not critical, it is preferred that the focal length is less than the distance from the lens of 
the transducer 1 0 to any solid material or to the opposing wall of a test cell, pipe or tank. In other words, the focal region 
5 should be within the fluid of interest. A transducer 1 0 having a focal length of about 1 00 mm to about 200 mm is suitable 
for the present invention as long as no wall, pipe or other process equipment structure intervenes between the trans- 
ducer 10 and the focal region. 

[0040] The selection of a transducer 1 0 involves consideration of its capability to introduce into the fluid an acoustic 
signal with respect to the characteristics of duration, power and frequency. Such signal can be in the form of a pulse 

10 that is transmitted to the transducer 10 as a high-voltage spike of short duration and typically repeated many times per 
second, or as a tone burst. Although the use of pulse signals is preferred, tone burst signals are useful in circumstances 
where the use of high voltage is undesirable and where a signal of narrower frequency band is preferred. The voltage 
of the input signal pulse can be from about 100 to about 600 volts. An input signal pulse voltage of from about 200 to 
about 500 volts is preferred and from about 250 to about 500 volts is more preferred. 

15 [0041 ] The input signal is transmitted to the transducer at a pulse repetition frequency of from about 5 to about 250 
Hertz (Hz). However, a frequency of from about 20 to about 200 Hz is preferred and a frequency of from about 100 to 
about 200 Hz is more preferred. 

[0042] The duration of the input signal can be from about 5 to about 80 nanoseconds (ns). However, it is preferred 
that the duration of the input signal be from about 10 to about 60 ns and it is more preferred that the duration be from 
20 about 20 to about 40 ns. 

[0043] The detected power of the input signal can be amplified so that it is from about 5 to about 50 decibels (dB). 
However, it is preferred that it be from about 10 to about 40 dB and more preferred that it be from about 20 to about 40 
dB. 

[0044] By way of example, a series of input pulses, each of 500 volts, having a duration of 25 nanoseconds and a 
25 power of 33 dB and repeated at a frequency of 200 Hz, and with a receiver gain of 33 dB, is suitable for some applica- 
tions of this invention. 

[0045] Typical operation of an acoustic transducer and ultrasonic systems similar to those suitable for use in the 
present Invention is described, for example, by Urick, R. J., J. AppL Rhys, 18, 983 - 987 (1947); McClements, D. J., et 
ai, J. Phys. D: Appl. Phys, 22, 38 - 47 (1989); Holmes, A. K., etal., J. Coll. Int. Sci, 156, 261 - 268 (1993); McClements, 
30 D. J.. Adv. ColL Int Sci, 37, 33 - 72 (1 991); McClements, D. J., J. Acoust Soc. Am, 91, 849 - 853 (1992); Pinfield, V. J,. 
et ai, J. Coll. Int. Sci, 166, 363 - 374 (1 994) and McClements, D. J., The use of ultrasonics lor characterizing fats and 
emulsions, Ph.D. Thesis, Food Science Department. University of Leeds, UK (1988). 

[0046] The transducer should be capable of resisting temperature to 200°C, and preferably to 300''C. and most 
preferably to 500°C. The transducer should also be capable of resisting pressure to about 200 kPa, and preferably to 
35 about 2.50 MPa and most preferably to about 5 MPa. Moreover, the transducer should preferably be capable of resisting 
chemical corrosion and physical erosion by the fluids in which it is to be used. 

[0047] Suitable transducers may be commercially obtained or may be constructed. One type of commercially avail- 
able acoustic transducer suitable for use in the present system is a SLIM 10 - 10/SF30 mm/834/02, 10 MHz transducer 
as supplied by Sonatest Ltd., Milton Keynes, U.K. 

40 [0048] The acoustic transducer 10 may also be used as a sensor to sense acoustic energy scattered by particles 
in the fluid. Alternately, a separate sensor 20 may be used. If a separate sensor 20 is used, it may be placed anywhere 
in the fluid in relationship to the probe 10, that is insensitive to the directly transmitted beam, but is preferably close 
enough to receive energy scattered from the focal region. The sensor, whether the transducer 10, or a separate sensor 
20, converts that portion of the scattered acoustic energy encountering the piezoelectric element into an electrical sig- 

45 nal. This signal con-esponds to the amplitude of the portion of scattered acoustic energy encountering the piezoelectric 
element as a function of time and Is termed the "output signal". Operation of an acoustic energy sensor, in general, is 
described by the references given above in the section describing the transducer 10. 

[0049] It is preferred that the sensor of the subject system has sensitivity suitable for sensing back-scattered energy 
having frequencies up to 20 MHz, but a sensor which can sense back-scattered energy having frequencies up to 100 
50 MHz is more preferred and a sensor whteh can sense back-scattered energy having frequencies up to 200 MHz is most 
prefenred. 

[0050] A pulser/receiver, or pulser 30 provides the necessary input signal to drive the transducer 10. An amplifier 
40 amplifies the sensed scattered signal before it is transmitted to the oscilloscope 50. The pulser 30 and amplifier 40 
may be separate components or may be combined in a single component. For example, a pulser/receiver that provides 
55 both pulser 30 and amplifier 40 functions In a combined component is a UTEX 320; available from Scientific Instru- 
ments Inc., Mississauqua, Ontario, Canada. It is preferred that the pulser/receiver is supplied with a software package 
such as, for example, Ultrasonic Instrument Control V3.04. The pulser/receiver is capable of providing an input signal 
having a frequency range of 1 MHz to 150 MHz; a pulse voltage of from 100 V to 500 V; a pulse width of from 5 ns to 
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80 ns; a pulse repetition frequency (PRF) of about 200 Hz to 20KHz; and a gain of from 0 to 63 dB 
[0051] The oscilloscope 50 of the present system should be capable of gating the scattered acoustic signal by lim- 
iting the output signal to the portion of acoustic energy scattered by particles in the focal region. The oscilloscope is also 
preferably capable of carrying out an analog-to-digital conversion and, more preferably, of transfomning the output sig- 
5 nal from an amplitude vs. time domain to a magnitude vs. frequency domain as described below. The gating function 
limits the output signal to the sensed energy that was scattered by particles within the focal region. Thus, the gating step 
deletes that part of the signal caused by any reflection due to an opposing pipe wall and all other portions of the signal 
except for the scatter due to particles in the focal region. 

[0052] An oscilloscope that is suitable for use in the present system is, for example, a LeCroy 931 OA, dual channel 
10 400 MHz oscilloscope, with a sampling rate up to 1 00 MS/s. Such oscilloscopes should be complete with software suit- 
able for waveform processing such as, for example, Wave form processing packages WP01 and WP02 incorporating 
FFT averaging function, available from LeCroy. 

[0053] Data from the oscilloscope is transferred to a computer 60 by using, for example, a National Instruments 
IEEE Plug and Play adapter. As an alternative, the oscilloscope may have an integral computer that is capable of car- 

15 rying out all of the necessary computational functions that will be described below. 

[0054] A computer 60 is preferably one that is capable of receiving the signal from the oscilloscope 50 and storing 
the waveforms for future reference and also of comparing the waveforms against reference scans and other standards. 
It is also preferred that the computer be capable of performing the derivation of the size, concentration and size distri- 
bution of particles from an acoustic speckle signal of the particles in a fluid; of determining the degree of asphaltene 

20 agglomeration in a hydrocarbon liquid and to initiate any desired alarm or control action. It is more preferred that the 
computer is programmed with one or more algorithms that are capable of relating the magnitude and/or duration of 
acoustic speckle signals to the size, concentration, and/or size distribution of the particles in the fluid. The algorithms 
can Include a diffusion model that is capable of relating the size of particles that participate in a speckle event to the 
duration of an acoustic speckle signal of the speckle event. As an alternative, the algorithm can be one that includes an 

25 autocon-elation model that is capable of relating the size of particles that participate in a speckle event to the duration 
of an acoustic speckle signal of the speckle event. 

[0055] While the type and computational speed of the computer are not critical, it Is prefen'ed that a personal com- 
puter having at least a Pentium® processor, or its equivalent, and having spreadsheet software such as, for example, 
Microsoft Excel®, be used in the present system. 

30 [0056] As briefly described above, the present invention provides a novel method for determining the size, concen- 
tration and size distribution of particles in a fluid. The method will now be described In more detail. 
[0057] The steps of the method include obtaining an acoustic speckle signal of the particles in the fluid and deriving 
the size, concentration and^orsize distribution of the particles from the acoustic speckle signal. One embodiment of the 
method for obtaining an acoustic speckle signal of the particles in a fluid is by introducing an Input signal of acoustic 

35 energy into the fluid; sensing a portion of scattered acoustic energy resulting from interaction of the particles with the 
input signal; producing an output signal corresponding to amplitude of the portion of scattered acoustic energy as a 
function of time; and transforming the output signal from a time domain to a frequency domain. A peak on the magni- 
tude versus frequency plot, as shown in Figure 5, can be identified as characteristic of the acoustic speckle. If the trans- 
formed output signal is averaged over a period of time that includes two or more speckle events, the coherent scattered 

40 acoustic energy will sum, the incoherent scattered acoustic energy will cancel out and a clearer acoustic speckle signal 
peak can be identified as shown by the less noisy line in Figure 5. 

[0058] The output signal that the oscilloscope receives is in an amplitude versus time domain. If it Is desired, the 
output signal can be limited to a certain frequency range before it is transmitted to the oscilloscope in order to delete 
high frequency signals that provide no Information about the particles In the fluid. In some situations, for example, the 

45 output signal can be limited to a frequency range of from 1 to 30 MHz, or even from 1 to 20 MHz. 

[0059] An example of an amplitude versus time signal on an oscilloscope for a sample of crude oil is shown in Rg- 
ure 6(a). In order to enhance the detection of acoustic speckle events, the output signal can be limited to acoustic 
energy scattered by particles in the focal region. Figure 6(b) shows another oscilloscope plot of an amplitude versus 
time output signal for crude oil that has been gated, or limited, to the signal from particles In the focal region. The 

50 selected part of the signal from the focal region is shown as a bolded section of the signal plot. 

[0060] The output signal from the focal region is then transfomned from a time domain to a frequency domain. This 
transformation is preferably carried out by the oscilloscope by using a Fast Fourier Transformation (FFT) algorithm to 
transform the gated signal from time to frequency domain. Computation of the FFT is described in detail at pages C-3 
to C-5 of the Operating Manual for the LeCroy Oscilloscope, LeCroy Corporation, Chestnut Hill, NY. Where a semlquan- 

55 tltative and simple assessment of the total backscattering is required, the power spectrum can be used. This provides 
a measure of the total backscattered intensity detected in the focal region. 

[0061] Figure 6(c) shows the oscilloscope plot of the FFT of the output signal of Figure 6(b). Figure 6(d) shows a 
plot of the FFT of the gated part of the transformed signal averaged over 900 oscilloscope sweeps. The time-averaged 
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acoustic speckle signal is identified on the plot of Figure 6(d) as the peak at a frequency of about 1 0 MHz - con^espond- 
ing to the transducer output maximum. 

[0062] As described previously, it is believed that the acoustic speckle signal for each individual speckle event has 
a certain duration that corresponds to the average size of the particles that produced the speckle event Once the dura- 
5 tion of an individual acoustic speckle event Is measured, the duration can be used to derive the size of the particles that 
produced the speckle event. One method for estimating particle size from the duration of the acoustic speckle signal is 
by the use of an algorithm derived on the basis of a simple diffusion model, such as: 

to = (Ay8)^/2D Equation 1 

70 

where: 

to is the duration of the acoustic speckle and Is the estimated time for dephasing of a acoustic mirror, in seconds. 
For this estimate the dephasing time is computed to be the time necessary for two particles to move 1/8 of a wave- 
rs length in opposite directions; 

X is the acoustic wavelength in meters; and 
D is the diffusion coefficient, in mV\ and Is 

D = kT/6;iTia Equation 2 

20 

where: 

k is the Boltzman constant (JK*^), (1 .38x1 0'^^); 
T is the absolute temperature in degrees Keh^in; 
71 is 3.14159; 
25 r| is the viscosity of the fluid in Pa • s; and 
a is the particle radius in meters; 

Thus, by combining the two equations and solving for the particle radius (a), the particle size can be calculated 
for each measurement of acoustic speckle duration (to). 
30 [0063] An alternative, and more accurate, diffusion model that can also be used for deriving particle size from a 
measurement of the duration of the speckle event Is: 

a = (t^)(4kT/67CTi)(8/Xe/fc(1))^ Equation 3 

35 where: 

a is the particle radius in meters; 

Is the duration of an acoustte speckle event In seconds; and other parameters are as defined above. 

40 [0064] As an alternative to deriving particle size from the duration of an acoustic speckle by using a diffusion model, 
the particle size can also be derived by using an autocon^elatlon model of the type described by Berry etaL, in Physical 
Chemistry, p.1106, Wiley Publ. New York (1980); and by Wellstead, RE., l\/1ethods and Applications of Digital Spectral 
Techniques, Technical Report No. 0008/83, UMIST, Solartron Instruments, Farnborough. U.K. 
[0065] Since the novel method measures the average size of particles that produce an individual speckle event, it 

45 is believed that when the duration of two or more speckle events are measured, the average partwie size that is calcu- 
lated for the different speckle events will vary within the range of the actual range of sizes of particles in the fluid. It fol- 
lows that the higher the concentration of particles of a given size, the more frequently those particles will produce 
speckle events. Thus, the concentration of particles of a given size can be obtained by measuring the frequency of 
occurrence of speckle events produced by particles having that given size. If desirable, one of ordinary skill would easily 

50 be able to obtain a correlation between particle size and frequency of occurrence of speckle events in a fluid having a 
given viscosity by carrying out a routine test in the same fluid with particles of known size at several concentration lev- 
els. 

[0066] The novel method is capable of making acoustic speckle measurements very rapidly and many speckle 
events can be recorded and measured within a short period of time. If a sufficient number of individual measurements 
55 are taken, the concentration of particles having average sizes that span the size range of all particles in the fluid can be 
obtained. This can be done by combining the particle concentration measurements for all of the average particle sizes 
measured. One method for making this combination, for example, is by defining the total range of particle sizes in the 
fluid as being the range between the smallest and the largest average particle size calculated in the test; dividing the 
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total range of particle sizes into a discrete number of smaller size ranges; and adding together the particle concentra- 
tions that were detemiined for particles having an average size within each smaller size range. For example, if the 
smallest and largest average particle sizes that were reported in the test were 1 micron and 1 00 microns, respectively, 
the total range of particle sizes in the fluid would be taken to be from 1 to 1 00 microns. This range could then be divided 

5 into almost any number of smaller size ranges - 20 for example. Each of the 20 smaller size ranges would, in this exam- 
ple, then include all particles of a size within a range of 100/20 = 5 microns. Thus, the discrete sub-ranges would be 0 
- 5 microns, 6-10 microns, and so on, up to 95-100 microns. An example of how such a plot of frequency of occur- 
rence versus average particle size might appear for this example is shown in Figure 7. The concentration of particles In 
the fluid could then be derived by adding the concentration of particles in each of the 20 discrete smaller size ranges. 

10 [0067] In fact, by using a similar concept, the novel method makes rt possible to derive the particle size distribution 
in the fluid from a frequency of occurrence versus average particle size distribution. Although only two or more of the 
particle size measurements are required for such derivation, it is preferred that many individual measurements be used. 
By way of example, if 2000 individual speckle events are measured over a period of 0.01 second, the 2000 speckle 
events can be categorized into a number of discrete average particle sizes, each discrete size having a nan-ow size 

IS range. The frequency of occurrence of speckle events In each discrete size range Is then essentially a measure of the 
size distribution of the particles in the fluid. In this case, the distribution measurement would have been provided in the 
short span of only about 0,01 seconds. This rapidity of measurement provided by the present invention would permit 
measurement and reporting of particle size distribution on a real-time basis. 

[0068] To obtain an accurate particle size distribution, a sufficient number of particle size calculations are obtained 
20 for the calculation of a particle size distribution for the particles in the fluid. If the plot of Figure 7, for example, is based 
on a sufficient number of individual particle size calculations, the curve itself would be indicative of the particle size dis- 
tribution. It is preferred that the size distribution be derived on the basis of at least about 1 00 individual particle size cal- 
culations and at least about 500 individual calculations are more preferred, at least about 1000 Individual calculations 
are even more preferred, and at least about 2000 calculations or more are most preferred. It is also necessary to pre- 
ss sume a particular fomn for the particle size distribution, such as, for example, a log-normal distribution. 

[0069] An alternative method for the calculation of the concentration of particles in a fluid comprises determining 
the size of the particles in the fluid from the duration of speckles as described previously - for example, by calculating 
the average radius (a) ~ and then detemiining the total cross-section area of particles in the focal region (Afr). It Is 
believed that the maximum magnitude (maximum peak height; H) of each individual speckle signal is directly propor- 
30 tional to the total cross-section area of the particles that produced the speckle event. Thus, by assuming that each par- 
ticle has a certain geometry (namely spherical), the total cross-section area of the particles that produced the Individual 
speckle event can be con-elated with the maximum magnitude of the individual speckle signal as: 

Afr = c(H) Equation 4 

35 

where: 

c is a proportionality constant; and 

H is proportional to ttie maximum peak height, or magnitude, of the individual speckle signal. 

40 

[0070] A representative volume can be defined as being the volume of fluid within the focal region (V^). If It is 
assumed that the particles that are within this volume are the particles that produce the acoustic speckle events that 
are being measured, tiien the concentration of particles in tiie focal volume can be derived from the volume of the focal 
region and the parttele sizes and total-cross section area of particles measured as: 

45 

n/Vjr = A|/jca^Vfr Equation 5 

where: 

50 nA/jr is the number of particles in the focal volume; and other parameters are as defined previously, 

[0071] If the concentration of particles in the focal region is measured for two or more speckle events and the values 
of such concentrations are averaged, tiie concentration of particles in the fluid can be calculated. As described previ- 
ously, it is preferred that measurements be taken for many speckle events in order to increase the accuracy of the cal- 
55 culation, 

[0072] For purposes of Increasing speed and simplicity of the method, the step of calculating particle concentration 
from speckle peak height can be left out and a value that is proportional to particle concentration can be cabulated 
directly from average particle size and speckle peak height data. 
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[0073] Once the concentration of particles of a given average size has been calculated by either nnethod, the size 
distribution of particles in the fluid can be calculated as described previously. 

[0074] The apparatus and method of the present invention can be used to measure the size, concentration and size 
distribution of particles in any fluid, but the novel apparatus and method are particularly useful for use where the fluid is 

5 a liquid and even more particularly when the liquid is an optically opaque liquid. The novel apparatus and method can 
be used when the liquid is corrosive, hot or under pressure. A particularly useful application of the subject method and 
apparatus is for the measurement of particles in crude oil, or other opaque hydrocarbon liquids. 
[0075] One embodiment of the novel method has been found to be particularly useful for determining the degree of 
agglomeration of small particles, such as asphaltene particles, in liquids. This embodiment is especially useful when 

10 the liquid is a hydrocarbon liquid. This method comprises obtaining, for two or more individual speckle events, an indi- 
vidual acoustic speckle signal from the hydrocarbon liquid and detennining the degree of asphaltene agglomeration 
from the acoustic speckle signal. 

[0076] In particular, an input signal of acoustic energy is introduced into the fluid, which input signal is focused on 
a focal region in the liquid. An output signal of the acoustic energy that is scattered by the particles in the focal region 
IS in response to their interaction with the input signal is sensed as amplitude as a function of time. This output signal is 
then transformed from a time domain to a frequency domain. As described above, the transformation is preferably car- 
ried out by an the oscilloscope by the use of a Fast Fourier Transform algorithm. An acoustic speckle signal from the 
transformed output signal in the frequency domain is then Identified. 

[0077] It should be emphasized that the common practice of averaging the instantaneous input signal will not work 
20 for the subject method, because the total of speckle events sums to zero and ali information that can potentially be 
derived from the speckle is lost. It is believed that this occurs because each speckle event occurs at a slightly different 
position in space. Thus, averaging of the bacl<scattered signal can only be earned out after the instantaneous signal 
amplitude versus time plot has been transformed into the frequency domain. 

[0078] Each of these steps is then repeated, but in a control liquid that contains no asphaltene particles. The trans- 
25 formed signal from the control liquid acts as a baseline standard. A plot of a transformed output signal from a crude oil 
sample plotted concurrently with a transfonned signal from n-heptane as a control fluid is shown in Figure 8. The 
degree of asphaltene agglomeration can be determined by comparing the acoustic speckle signal from the hydrocarbon 
liquid with the acoustic speckle signal from the control liquid. 

[0079] One method by which the signal from the hydrocarbon liquid and the sample from the control liquid can be 
30 compared is by subtracting the magnitude of the acoustic speckle signal of the hydrocariaon liquid at a selected fre- 
quency from the magnitude of the acoustic speckle signal of the control liquid at the same frequency to obtain a differ- 
ence. The maximum value of this difference between the two curves is herein termed the "speckle heighf. An 
illustration of the calculation of the speckle height from the data of Figure 8 is shown in Figure 9, where maximum 
speckle height is seen to occur at a frequency of about 1 1 MHz. 
35 [0080] The method and apparatus of the present invention can be used when the test fluid is static, such as in a 
non-flow sample cell, but it can also be used in an agitated tank or process line In which the fluid is flowing. The only 
limitation for the use of the method in a flow situation is that the partfcles constituting the speckle mirror remain illumi- 
nated by the acoustic beam throughout the speckle event. Thus, flow along the acoustic beam axis will remain longer 
in the beam than flow transverse to the beam. In practice, the longest event that can be detected will be a function of 
40 beam width, flow rate, the size of the focal region and the relationship between the direction of fluid flow and the acous- 
tic beam direction. Transducers that are installed in an elbow of a flow conduit and which have their beam aligned into 
the incoming or outgoing fluid, for example, would be preferred for use in fluid flow situations. 
[0081] The present method and apparatus have the advantage of being useable with only one transducer when it 
is important to minimize cost or number of penetrations into the fluid. Only a single type of signal is required - an accus- 
es tic energy signal - rather than requiring other, additional signals and there is no need for the addition of titrants, such 
as heptane or hexane, or contaminating chemicals to the subject fluid. I\4easurements by the present method and appa- 
ratus provide a measurement of parttele size, concentration and/or size distribution that are available very quickly and 
even on a real-time basis. 

[0082] The following examples describe the prefen-ed embodiments of the invention. Other embodiments within the 
so scope of the claims herein will be apparent to one skilled in the art from consideration of the specification or practfce of 
the Invention as disclosed herein. It Is intended that the specification together with the examples, be considered exem- 
plary only, with the scope and spirit of the invention being Indicated by the claims which follow the examples. 

EXAMPLE 1 

55 

[0083] This example shows the observation of acoustic speckle in crude oil. 

[0084] A SLIM 10-10/SF30mm/834/02 transducer, available from Sonatest, UK; Dickens Road, Milton Keynes, 
MK12 5QQ, UK; was connected with a pulser/receiver and a LeCroy 931 OA oscilloscope as shown in Figure 3. The 
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transducer also acted as the sensor for backscattered acoustic energy. The input signal settings were 500 volts, 25 ns, 
200 Hz and using a power receiver having a gain of 33 dB. The oscilloscope settings were set to an Impedance of DC60 
Ohms, with the screen plot axes being set to 0.5 nnicroseconds per division for the X-axis (time) and from 5 mV to 20 
mV per division (amplitude) for the Y-axis. 

5 [0085] The tip of the transducer was immersed into a sample of crude oil at 23*'C in a test cell and the input signal 
and sensing of the scattering was started. The delay time required for the input signal to move to the focal region and 
the backscattered acoustic energy to move back to the transducer was 41 .5 microseconds + 0.5 microseconds/division 
X 5 divisions. 

[0086] The backscattered signal received by the transducer was gated to the focal region and transfonned by the 
70 oscilloscope by using a Fast Fourier Transfonn from a time domain to a frequency domain. The FFT transformed signal 
was averaged over 435 sweeps and the magnitude of the resulting averaged transformed signal was plotted versus fre- 
quency. A peak was seen which had a maximum magnitude at about 10-12 MHz. This peak was identified as being 
the "speckle peak". It is believed that the peak is the Image in the frequency domain of the acoustic transducer output. 

75 EXAMPLE 2 

[0087] This example illustrates the calculation of particle size as a function of the duration of an acoustic speckle 
by the use of two diffusion models. The first model Is a simple calculation, but provides a less accurate value than the 
second model. The physical constants for asphaltene particles in crude oil were estimated and used for the computa- 
20 tion for both models. The second model is used to show the dependence of the calculated particle size on the viscosity 
of the fluid. 
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Definition of variables 


PARAMETER 


UNITS 


SYMBOL 


VALUE 


FORMULA OR COMMENTS 


Boltzman constant 




k 


1.38x10" 




Tenftperature 


""C 




25 


Ambient 


Absolute temperature 


K 


T 


298.15 




Vrscosity 


Pas 


n 


7.00x10-^ 


Figure for low viscosity fraction 


Particle Radius 


m 


a 


1.00x10"* 


Figure for asphaltenes 


Particle Dianneter 


m 


d 


2.00x10-* 


=2a 


Particle density 




P 


1.00x10^ 


Guess 


Particle Mass 


kg 


m 


4.19x10'^' 














Velocity of sound 


ms*^ 


V 


1200 


Approximate figure for crude oil 


Frequency of sound 


Hz 


f 


1.00x10^ 


Sonatest transducer 


Wavelength 


m 


X 


1.2x10*^ 





Computed values of the duration of a speckle event for a particle with radius of 10** m 

Mobility coefficient' 



skg' 
m s 



Diffusion coefTtcient 



Tlnie interval** 



O 



2.50x10^ 
3.12x10"^' 



36.1 



=5fl/d 
kT 



r {\JBfl2D 



A more accurate version of a diffusion model can be formulated as shown below: 





-1.05x10'^^ 




1.05x10-^^ 


D(t) 


3.12x10*'^ 



= -mn + exp(-to/m|j) 

=f77fJ 

=D(1.(tcoa/m^)) 



Time correction s 
Einstein condition^ 
(fo»rn//) 

Time corrected diffusion m s 
coefficient*^ 

The above three lines are checks on the validity of the Einstein condition and include on the 
final line a corrected version of the diffusion coefficient in the case that the Einstein condition 
fails. This is not needed in the example calculation given here but can be necessary under 
some conditions. 



<r> 



2.37x10"^ 



,0.5 



0.446 



' ((X/8)erfc(1))*/4D 



Distance moved in time m 
inteivar 

Period s T(0 1.0x10'' 

Time for a step change s 
in concentration to 
reduce to 50% over a 
distance of XJS 

The above is a more accurate caiculatk)n of the time taken for the particle mirror to break up 
due to diffusion. It agrees very well with experimental obsen^attons and with what is known 
about the size of asphaltene particles. 

Dickinson, E. and McCIements, DJ., Advances in Food Colloids*\ Blackie, 1995, 
pi 26, equation 4.2S, estimate for particle in water. Here it is assumed that the 
dififusion of asphaltenes is similar to that of particles in water. 
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• This is the condition for the de-phasing of the mirror. It is derived from the Einstein 
equation for the average displacement of a given particle in a one-dimensional 

5 

projection (Einstein, 1905), and see reference above, page 125, equation 4.23. The 
condition for de-phasing of the mirror is that particles need only move 1/8 of a 
wavelength. A movement of one-half of a wavelength is necessary for complete phase 
cancellation. However, at least two particles must be involved and they can move in 
opposite directions. So a movement of one-eighth of a wavelength on the part of two 
particles moving in opposite directions, along the transducer axis, will result in a 
substantial reduction of the phase of the mirror. 

^' Equation 4.25 in Dickinson, E. and McClements, D.J., Advances in Food 
20 Colloids'\ Blackie (1995). 

This is a different way of writing Equation 4.25 in Dickinson and McClements, Id. 
This is computed from the Einstein equation for motion in 3 dimensions. <r^> = 

The viscosity of water is assumed to be 0.008 Pa-s; and the viscosity of sunflower oil 
is 0.05 Pa s. 

30 

[0088] TTie table below is in three parts. The first section of the table shows the effect of variation in the viscosity of 
the fluid on the particle size calculation. The second section of the table shows the dependence of the duration of an 
35 acoustic speckle on particle size at a single value of viscosity. The third part of the table shows the calculation of particle 
size as a function of the duration of an acoustic speckle by a diffusion model. 
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5 



25 



Dependence of t^irror ^ 


V 


^mirror 


D 


a 


0.0001 


0.006376 


2.18279x10-1° 




1.00x10'8 




0.001 


0.063762 


2.18279x10"^^ 




1.00x10*® 




0.01 


0.637622 


2.18279x10-12 




1.00x10-8 




0.1 


6.376218 


2.18279x10 ''^ 




1.00x10-8 




Dependence of tfnffror ^ 


a 


^mirror 


D 




1.00x10'^ 


0.446335 


3.12x10-12 




7.00x10'^ 




1.00x10-^ 


4.463352 


3.12x10-1^ 




7.00x10'^ 




1.00x10"® 


44.63352 


3.12x10-""* 




7.00x10'^ 




1.00x10-5 


446.3352 


3.12x10'i5 




7.00x10-^ 




I.OOxlO-"^ 


4463.352 


3.12X10-1® 




7.00x10'^ 




1.00x10'^ 


44633.52 


3.12x10-1^ 




7.00x10-3 




Dependence of a on f,ni^n>r 




a ^meters) 




1.00x10-"^ 


2.24x10-12 








1.00x10-^ 


2.24x10-11 




(4kT/67iTi)((8/Xerfc(1))^) 


1.00x10-2 


2.24x10-1° 








LOOxlO"" 


2.24x10"® 








1 




2.24x10-8 








10 
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[0089] All references cited in this specification, including without limitation all journal articles, brochures, manuals, 
40 periodicals, texts, manuscripts, website publications, and any and all other publications, are hereby incorporated by ref- 
erence. The discussion of the references herein is intended merely to summarize the assertions made by their authors 
and no admission is made that any reference constitutes prior art Applicants reserve the right to challenge the accu- 
racy and pertinence of the cited references. 

[0090] In view of the above, it will be seen that the several advantages of the invention are achieved and other 
45 advantageous results are obtained. 

[0091] As various changes could be made in the above methods and compositions without departing from the 
scope of the invention, it is intended that all matter contained in the above description and shown in the accompanying 
drawings shall be interpreted as illustrative and not in a limiting sense, distribution of small particles in optically opaque 
liquids. 

50 

Claims 

1 . A method for determining the size of particles in a fluid comprising the steps of obtaining an acoustic speckle signal 
of the particles in the fluid and deriving the size of the particles in the fluid from the acoustic speckle signal. 

55 

2. A method for measuring the concentration of particles in a fluid comprising the steps of obtaining an acoustic 
speckle signal of the particles in the fluid and deriving the concentration of particles in the fluid from the acoustic 
speckle signal. 
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3. The method as set forth in claim 2 wherein the acoustic speckle signal is obtained by introducing an input signal of 
acoustic energy into the fluid, the input signal being focused on a focal region in the fluid and some of the input sig- 
nal being scattered by particles in the focal region; sensing a portion of scattered acoustic energy resulting from 
interaction of the particles with the input signal; producing an output signal limited to the portion of acoustic energy 
scattered by particles in the focal region and corresponding to amplitude of the portion of scattered acoustic energy 
as a function of time; transforming the output signal from a time domain to a frequency domain; and identifying the 
acoustic speckle signal from the output signal transformed to the frequency domain; and wherein the concentration 
of particles in the fluid is derived by determining the size of the particles in the fluid; determining the total cross- 
section area of particles In the focal region; detennining the volume of the focal region; and calculating the concen- 
tration of the particles in the fluid from the size of the particles, the total cross-section area of the particles in the 
focal region and the volume of the focal region. 

4. The method as set forth in claim 3 wherein the fluid contains particles that produce individual speckle events and 
wherein the concentration of the particles in the fluid is derived from the acoustic speckle signal by, for each of two 
or more individual speckle events, identifying an individual acoustic speckle signal con-esponding to each individual 
speckle event; measuring the duration of the individual acoustic speckle signal; calculating the size of the particles 
that participate in the individual speckle event from the duration of such individual acoustic speckle signal; deter- 
mining the total cross-section area of the particles that produced the individual speckle event by measuring the 
maximum magnitude of each individual speckle signal and correlating the total cross-section area of the particles 
that produced the individual speckle event with the maximum magnitude of the individual speckle signal; and deter- 
mining the concentration of particles in the fluid from the volume of the focal region and the particle sizes and total- 
cross section area of particles measured for the two or more speckle events. 

5. A method for detenriining the size distribution of particles in a fluid comprising the steps of obtaining an acoustic 
speckle signal of the particles in the fluid and deriving the size distribution of the particles from the acoustic speckle 
signal. 

6. The method as set forth in any of claim 1 , 2 or 5 wherein the acoustic speckle signal is obtained by introducing an 
input signal of acoustic energy into the fluid; sensing a portion of scattered acoustic energy resulting from interac- 
tion of the particles with the input signal; producing an output signal corresponding to amplitude of the portion of 
scattered acoustic energy as a function of time; transforming the output signal from a time domain to a frequency 
domain; and identifying the acoustk: speckle signal from the output signal transfomied to the frequency domain. 

7. The method as set forth in claim 6 wherein the input signal is focused on a focal region in the fluid and some of the 
input signal is scattered by particles in the focal region, and wherein the output signal is limited to the portion of 
acoustic energy scattered by particles in the focal region. 

8. The method as set forth In claim 7 wherein the fluid contains particles that produce individual speckle events and 
wherein the size distribution of the particles in the fluid is derived by calculating, for two or more individual acoustic 
speckle events, the size of the particles that produced each Individual speckle event; measuring the frequency of 
occurrence of speckle events that were produced by particles having the same size; and deriving the particle size 
distribution from the frequency of occun-ence of speckle events as a function of parttele size. 

9. The method as set forth in claim 8 wherein a sufficient number of particle size cateulations are obtained for the cal- 
culation of a particle size distribution for the particles in the fluid and wherein the size distribution of particles in the 
fluid is detennined as the frequency of occun-ence of acoustic speckle events produced by particles having a cer- 
tain size. 

10. The method as set forth in claim 9, comprising at least about 1 GO particle size calculations, preferably at least about 
500 particle size calculations, more preferably at least about 1 000 particle size cateulations. and most preferably at 
least about 2000 particle size cateulations. 

11. The method as set forth in claim 7 wherein the fluid contains particles that produce individual speckle events and 
wherein the size of the particles in the fluid is detennined from the acoustic speckle signal by, for each of two or 
more individual speckle events, Identifying an individual acoustic speckle signal corresponding to each individual 
speckle event; measuring the duration of the Individual acoustic speckle signal; deriving the size of the particles 
that produced the speckle event from the duration of the individual acoustic speckle signal; and determining the 
size of particles in the fluid from the particle sizes measured for the two or more speckle events. 
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12. The method as set forth in claim 1 1 wherein, for each individual speckle event, the size of the particles that pro- 
duced the speckle event is derived from the duration of that individual speckle event by using an autocon-elation 
model. 

13. The method as set forth in claim 11 wherein, for each individual speckle event, the size of the particles that pro- 
duced the speckle event is derived by con'elating the particle size with the maximum amplitude of such acoustic 
speckle signal. 

14. The method as set forth in claim 7 wherein the fluid contains particles that produce individual speckle events and 
wherein the particle concentration is derived by determining, for two or more Individual acoustic speckle events, the 
size of the particles that produced each individual speckle event; measuring the frequency of occun-ence of speckle 
events that were produced by particles having the same size; correlating the frequency of occurrence with partfcle 
concentration for each particle size determined: and combining the particle concentration for all particle sizes 
determined to arrive at the concentration of particles in the fluid. 

15. The method as set forth in claim 7 wherein the fluid contains particles that produce individual speckle events and 
wherein the size distribution of the particles in the fluid is detemiined from the acoustic speckle signal by, for each 
of two or more individual speckle events, identifying an individual acoustic speckle signal corresponding to the indi- 
vidual speckle event; measuring the duration of the individual acoustic speckle signal; deriving the size of the par- 
ticles that produced the individual speckle event from the duration of the individual acoustic speckle signal; 
determining the volume of the focal region; calculating the concentration of particles having the size of those that 
produced the individual speckle event from the particle size and the volume of the focal region; and determining the 
size distribution of particles in the fluid from the particle size and concentration as measured for the two or more 
speckle events.. 

16. The method as set forth in any of claim 4, 8, 11 . 14 or 15 wherein, for each individual speckle event, the size of the 
particles that produced that speckle event is derived from the duration of the Individual acoustic speckle signal for 
that individual speckle event by using a diffusion model. 

17. The method as set forth In claim 16 wherein the diffusion model is of the form: 

a=: (t^)(4kT/6jni)(8/Xerfc(1))^ 

where: 

a is the particle radius in meters; 

tm Is the duration of an acoustic speckle event In seconds 

X is the acoustic wavelength in meters; 

k is the Boltzman constant (JK'^), (1 .38x1 Q-^^); 

T is the absolute temperature in degrees Kelvin; 

K is 3.14159; and 

1] is the viscosity of the fluid In Pa • s. 

18. The method as set forth in any of claim 4, 8, 14, or 15 wherein, for each individual speckle event, the size of the 
particles that produced that speckle event is derived from the duration of the individual acoustic speckle signal for 
that individual speckle event by using an autocorrelation model. 

19. The method as set forth in any of claim 3, 8, 14 or 15 wherein, for each individual speckle event, the size of the 
particles that produced that speckle event is derived by measuring the amplitude of the individual acoustic speckle 
signal and correlating the particle size with the maximum amplitude of such individual acoustic speckle signal. 

20. The method as set forth in claim 5 wherein the fluid is a liquid, preferably a hydrocariaon, more preferably an opti- 
cally opaque hydrocarbon liquid, most preferably crude oil. 

21. The method as set forth in claim 20 wherein the input signal of acoustic energy is an Input signal of ultrasonic 
energy. 

22. The method as set forth in claim 21 wherein the particles are asphaltene particles. 
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23. The method as set forth in claim 1 or 2 wherein the fluid is crude oil and the particles are asphaltene particles. 

24. A method for determining the degree of particle agglomeration in a liquid containing particles that produce individ- 
ual speckle events, the method comprising obtaining, for two or more individual speckle events, an individual 

5 acoustic speckle signal from the liquid and determining the degree of agglomeration from the acoustic speckle sig- 
nal. 

25. The method as set forth in claim 24 wherein the particles are asphaltene partteles and the liquid is a hydrocarbon 
liquid. 

10 

26. The method as set forth in claim 25 wherein the method comprises the steps of introducing an input signal of 
acoustic energy into the fluid, which input signal is focused on a focal region in the liquid; sensing as an output sig- 
nal the amplitude as a function of time of the acoustic energy that is scattered by the particles in the focal region in 
response to their interaction with the input signal; transforming the output signal from a time domain to a frequency 

75 domain; identifying an acoustic speckle signal from the transfomied output signal in the frequency domain; repeat- 
ing each of the steps listed above, except in a control liquid that contains no asphaltene particles; and determining 
the degree of asphaltene agglomeration by comparing the acoustic speckle signal from the hydrocarbon liquid with 
the acoustic speckle signal from the control liquid. 

20 27. The method as set forth in claim 26 wherein the control liquid is n-heptane. 

28. The method as set forth in claim 25 wherein the magnitude of the acoustic speckle signal of the hydrocarbon liquid 
at a selected frequency is subtracted from the magnitude of the acoustic speckle signal of the control liquid at the 
same frequency to pbtain a difference. 

25 

29. The method as set forth in claim 28 wherein the maximum difference between the acoustic speckle signal of the 
hydrocarbon liquid and the acoustic speckle signal of the control liquid is identified as a speckle height. 

An apparatus for determining the size distribution of particles in a fluid comprising a signal generator capable of 
generating an electrical signal; a transducer that can translate the electrical signal from the signal generator into an 
acoustic input signal that is introduced into the fluid; a sensor that can sense acoustic energy that is scattered by 
the particles in the fluid in reponse to their interation with the acoustic input signal and provide an output electrical 
signal that is related to the sensed acoustic energy; an oscilloscope that is capable of acquiring the output signal 
in an amplitude versus time domain and transfomiing the signal to a magnitude versus frequency domain; and a 
computer that is programmed with an algorithm that is capable of identifying an acoustic speckle signal in the fre- 
quency domain and relating the magnitude and/or duration of the acoustic speckle signal to the size distribution of 
the particles in the fluid. 

31. The apparatus as set forth in claim 30 wherein the transducer has a focussing lens. 

40 

32. The apparatus as set forth in claim 30 wherein the algorithm comprises a diffusion model that Is capable of relating 
the size of particles that participate in a speckle event to the duration of an acoustic speckle signal of the speckle 
event. 

45 33. The apparatus as set forth in claim 30 wherein the algorithm comprises an autocorrelation model that is capable of 
relating the size of partbles that participate in a speckle event to the duration of an acoustic speckle signal of the 
speckle event. 
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PIGURE 1(b). 
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FIGURE 2(a). 
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FIGURE 2(b). 
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FIGURE 2(c). 
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FIGURE 6(a). 
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